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ABSTRACT 

Using the potential of two unprecedented missions, STEREO and RHESSI, 
we study three well observed fast CMEs that occurred close to the limb together 
with their associated high energy flare emissions in terms of RHESSI HXR spec- 
tra and flux evolution. From STEREO/EUVI and STEREO/COR1 data the full 
CME kinematics of the impulsive acceleration phase up to ~4 R Q is measured 
with a high time cadence of <2.5 min. For deriving CME velocity and acceler- 
ation we apply and test a new algorithm based on regularization methods. The 
CME maximum acceleration is achieved at heights h < 0.4 R , the peak velocity 
at h < 2.1 R Q (in one case as small as 0.5 R©). We find that the CME accelera- 
tion profile and the flare energy release as evidenced in the RHESSI hard X-ray 
flux evolve in a synchronized manner. These results support the "standard" 
flare/CME model which is characterized by a feed-back relationship between the 
large-scale CME acceleration process and the energy release in the associated 
flare. 

Subject headings: Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: hard 
X-rays 
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Introduction 



Solar flares and coronal mass ejections (CMEs) are the most violent phenomena in our 
solar system. Many aspects of the basic physics of these events are still not well understood. 
In the "standard" model it is envisaged that the erupting filament or CME stretches the 
coronal magnetic field lines to build up a vertical current sheet, where magnetic reconnection 
sets in to explosively release vast amounts of free magnetic e nergy, previou sly stored in the 



Forbes 



20001). The released 



corona in non-potential magnetic fields (for a review see e.g. 
energy goes into plasma heating, acceleration of particles to suprathermal velocities as well 
as into kinetic energy of the eruption. In this model, a close relation between the kinematics 
of the CME and the energy release in the associated flare is expected. 

A signific ant fraction o 



electrons (e.g. 



Hudson et al. 



the p rimary flare energy goes directly into acceleration of fast 
19921 ). As the accelerated electrons precipitate downward along 
the newly closed magnetic field lines to the lower lying denser atmospheric layers, they are 
collisionally stopped and heat and ionize the chromosphere and lower transition region. 
This can then be observed as enhanced UV and Ha radiation. If the beam flux is high 
enough, they will also emit detectable hard X-rays (HXRs) via nonthermal bremsstrahlung 
when the electrons scatter off ions of the ambient thermal plasma. Thus, HXR emission 
provides the most direct indicator of the evolution of the energy release in a flare (e.g. 



Fletcher fc Hudson 



20011). 



Several authors 



soft X-ray emission ( 


Zhang et al. 


2001 


Maricic et al. 


2007 




Vrsnak et al. 


2007) 



bund a correlation between the CME acceleration and the flare 



2001 



Zhang fc Dere 



2006 



Vrsnak et al. 


2004 




Temmer et al 


( 


2008) 



20071 ). In a recent case study by 
an almost synchronized behavior between CME acceleration and flare HXR emission was 
obtained for two on-disk events. Such results reveal that, for some associated flare-CME 



events, the reconnection process during the flare is closely related to the CME kinematical 
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evolution. The question remains if it is also possible to relate the flare HXR emission 
(count rate, spectral parameters) to the CME accele ration magnitu de. The Reuven Ramaty 



Lin et al. 



20021 ) delivers HXR spectra 



High-Energy Solar Spectroscopic Imager (RHESSI; 
and images at high temporal and spectral resolution, which enables us to study in detail 
the flare energy release process. 

To study the flare energy release in relation to the CME kinematics, the impulsive or 
main acceleration phase of a CME has to be cover ed. Since the impulsive ac c eleration phas e 



of a CME ta 



1999 



tes pl a ce at distances of R < 3 R^ (IMacQueen fc Fisher 



Vrsnak 



2001 



Zhang et al. 



2001 



Temmer et al 



1983 



St. Cyr et al 



20081 ) observations at low coronal 



heights and of high temporal cadence are required but are only limited available from white 
light coronagraphs. The Extreme Ultraviolet Imager instr uments aboard the twin spacecraft 



Kaiser et al. 



20081 ) mission have 



of the Solar Terrestrial Relations Observatory (STEREO; 
a large field of view and a high time resolution in the 17lA passband, perfectly suiting such 
studies. Together with COR1 observations, the inner coronagraph aboard STEREO, the 
impulsive acceleration phase of a CME is fully covered from its launch in the low corona up 
to 4.0 R . 

In this paper we study and analyze the CME dynamics and HXR emission of the 
associated flare for three well observed CME/flare events, using the potential of two 
unprecedented missions, STEREO and RHESSI. The events are selected to be located 
close to the limb in order to minimize the effect of projection in the CME kinematics. 
In addition, we apply a new algorithm (based on regularization methods) to derive the 
CME velocity and acceleration from the measured height-time data. A systematic test 
of the regularization method as well as a least-squares spline algorithm to a variety of 
synthetic CME acceleration profiles is performed, in order to evaluate the limitations and 
uncertainties of both methods. 
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Data and Methods 



The Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI; 



Howard et al 



STEREO-A (head) and ST 
Ultravi olet Imager (EUVI; 



20081 ) instrument package is part of each of the STEREO twin spacecraft, 
REO-B(ehind). It includes among others the Extreme 
Wuelser et al. 200^ and the white-light inner coronagraph 



COR1 ([Thompson et al. 



2003|). EUVI has a field of view (FoV) of 1.7 Rq which enables 
us to follow the erupting CME structure during its initiation and early propagation phase. 
Combining EUVI with COR1 observations, which have a FoV of 1.4-4.0 R Q , we can derive 
a complete velocity and acceleration profile for the early CME evolution. The partially 
overlapping FoVs between EUVI and COR1 allow us to check if the same features are 
followed in both instruments. 

Further advantages are the high time cadence of EUVI and COR1 observations which 
allow us to study the evolution of impulsive CME events low in the corona. EUVI observes 
in the EUV 17lA passband (Fex,xi: T~ 1 x 10 6 K) with a nominal cadence of 2.5 min; 
COR1 has a cadence of 10 min. Each day from 18:00 to 22:00 UT, EUVI and COR1 gather 
observations with increased cadence, namely 75 sec for EUVI and 5 min for COR1 in order 
to coordinate with Mauna Loa Solar Observatory as well as during observing campaigns. 
For the study of the CME kinematics, we combine EUVI 17lA filtergrams and total 
brightness images calculated from COR1 polarization sequence triplets. The level-0 data 
are properly reduced using SolarSoft routines, all images are background subtracted (for 
COR1 monthly background is subtracted from each polarization component) and rotated 
to solar north up. 



RHESSI performs imaging spectroscopy of solar flares in t he energy 



to 17 MeV with high temporal, spatial and spectral resolution (ILin et al. 



range from 3 keV 



2002|). The flares 



under study are fully covered by RHESSI observations during their peak phase which 
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enables us to study in detail their HXR spectra and flux evolution in relation to the 
kinematics of the associated CMEs. RHESSI spectra were integrated over 20 sec during the 
flare peak with a spectral resolution of 1 keV and fitted w ith a thermal plus non-thermal 



power-law component using OSPEX ( iSchwartz et al. 



2002) • From the power-law component 



we derive the amplitude of the HXR photon spectrum at 50 keV, the spectral index 7 of the 
photon spectrum (slope of the power law), and th e power in electrons above a cutoff energy 



197ll ). RHESSI images wer e reconstructed 



of 25 keV for thick-target emission (IBrown 
during the flare peak using the CLEAN algorithm flHurford et al. 
detectors 2 to 8 (except 5 and 7). 



2002j) including front 



3. Derivation and test of CME acceleration profiles 



For our study the reliable derivation of CME velocity and acceleration profiles from 
the measured height-time (HT) data is crucial. The central problem is to properly smooth 
the noisy data and to estimate the impact of measurement errors on the derived quantities, 
especially on the acceleration profile. This issue receive d a lot of attent ion for deriving 



CME kinematical curves from SOHO/LASCO data (see 



Wen et al. 



20071 ). Here, we apply a 



new type of technique to derive velocity and acceleration profiles fro m HT measurements 



Kontar et al. 



(120041 ) 



namely an enhanced regularization algorithm originally developed by 
to invert solar X-ray spectra measured by RHESSI. 

The main problem of time derivatives of data with measurement errors is that its finite 
difference estimate always leads to an amplification of the error. This ca n be expressed 



as the sum of discretization (finite time cadence) and propagation errors (IGroetsch 



1984h 



The former error is proportional to the time cadence of the measurements, while the latter 



i s inversely proportional to it. There 



(IHanke fc Scherzer 



200 lh . Following 



ore, the total deri v ative error always has a minimum 



Hanke fc Scherzer! (120011 ). and 



Kontar fc MacKinnon 
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( 120051 ) the regularized method searches for a model-independent velocity and acceleration 
estimate (regularized solution) with a minimum error, which produces smooth derivatives 
an d avoids additional err o rs typ ical of finite differences. In addition, the regularized method 



of iKontar fc MacKinnon! (j2005l ) provides the confidence interval for the first and second 
derivative, i.e. for CME velocity and acceleration profiles. The outcome of the regularization 
algorithm is controlled by two parameters. The number of time bins (which has a smoothing 
effect) and the regularization "tweak", an important parameter that regulates how reliable 
are the uncertainties of the input data (e.g. "tweak" of one means that the errors follow a 



normal distribution without systematics or correlation). 



J 



To test the reliability of the results we compare the CME acceleration profiled derived 
by the regularization algorithm with those derived from a least-squar es spline fit methoc 



which is an advanc ed method already used in previous CME studies (IMaricic et al. 



Vrsnak et al 



2004 ; 



20071 ). (Note that we did not compare less sophisticated methods, like direct 



derivation of the HT data or 3-point Lag rangian interpola t ion si nce these would give worse 



results.) Based on the formulas given in 



Gallagher et al. 



(120031 ). we generate synthetic 



kinematical curves of CME propagation representing different scenarios of CME evolution 
(impulsive and gradual acceleration profiles). 

We then reduce the sampling rate of the generated curve and add some noise. The 
noise consists of random noise as well as of manually shifting individual data points. With 
this we aim to take into account the decrease in the dynamic range of CME observations 
from the inner to the outer corona as well as inconsistencies of tracking features across 
different wavelength regimes (EUV and white light). From these "imperfect" data points 



deriving the acceleration profile, i.e. calculating the second derivative, highly intensifies 
noise on the data. The outcome of the acceleration is critical in order to reliably compare 
with HXR emission. 
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the acceleration profile is derived and compared to the acceleration profile from the original 
"perfect" synthetic curve using the regularization and the spline method. To test the 
stability of the two methods, we simulate HT data with different errors and sampling rates. 
For both methods the normalized residuals are derived, i.e. the differences between the 
measured HT values and the values resulting from the fitting procedure and the regularized 
solution, respectively, divided by the error. 

In addition to the regularization method we use the spline fit routine which computes 
least squares splines with equally spaced nodes to the HT data points (SPLFIT in IDL). 
From the spline fit curve, we numerically derive the second order derivative to determine 
the accelera tion. The smoothin g; effect of the spline fit is controlled by the number of 



nodes n (cf. 



Vrsnak et al. 



20071 ). To illustrate the way in which this affects the shape of 
the derived acceleration curve, we present the results of the spline fit for three different 
numbers of nodes (n — 1, n, n + 1). The nodes are chosen in such a way that the residuals 
for each of the curves would be equally small. 

We first test the reliability of the methods and derive the acceleration from the 
synthetical HT data without noise (Fig. [TJ left panels). For comparison, the right panels of 
Fig. [1] show the results derived from HT data with noise added. We note that the second 
part of the synthetic curve (decreasing height) is unrealistic for a CME evolution, but the 
profile was used to provide a challenge to both methods. As can be seen from the left 
panels, without noise the outcome from the regularization tool represents the acceleration 
curve better than that from the spline fit with respect to the timing of the acceleration peak 
and the acceleration duration, though the spline fit method reproduces the HT curve quite 
well. The peak amplitude of the CME acceleration is well represented by the regularization 
tool, whereas it is underestimated (by about 20%) by the spline fit method. In addition, the 
regularization tool reveals horizontal error bars (±1.0 min) to account for the uncertainty in 
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timing due to the reduced sampling rate. From the right panels of Fig. [I] it can be seen that 
the acceleration curve derived from the spline fit applied on noisy data has only slightly 
changed, whereas those from the regularization tool did since it considers the information 
provided by the error bars to yield an adequate error estimation which we do not get from 
the spline fit method. Strictly speaking, as soon as we add noise to the true HT data, the 
"true" velocity or acceleration cannot be reproduced. So the prime objective of any method 
is to provide the range of velocities/accelerations where the "true" solution should be. 

Figure [2] shows the outcome for a HT curve of an impulsively accelerated CME 
evolution with large noise added (average error of 0.33 R , i.e. larger than typically derived 
from real observations). We see that the acceleration peak is neither well represented by 
the regularization tool nor by the spline fits. However, the intrinsic errors in time which we 
get from the regularization tool (±2.0 min) include the true solution. The overall shape, 
i.e. the duration of acceleration, is well represented by both methods. 

Figure [3] (left) shows a CME evolution curve with small errors (±0.07 R ) and a 
sampling rate increasing from 2 to 10 minutes, simulating a cadence of measurement points 
comparable to the actual observations in the present study. The spline fit acceleration varies 
only little when applying different numbers of nodes (6,7,8) and matches well the peak but 
not the duration of acceleration. The outcome of the regularization tool for the most part 
matches the true solution within its error bars (uncertainty in the timing of the acceleration 
peak about ±3.0 minutes). The right panel of Fig. [3] shows the same kinematical curve, 
however, with a higher sampling rate of 1 and 2 min, respectively. The acceleration peak 
derived from the spline fit using different nodes shows a larger discrepancy. The timing 
of the acceleration peak calculated from the regularization tool is consistent with the true 
curve. However, its uncertainty is estimated to about ±4.0 min. To summarize, a sampling 
rate of ~2 min covering the impulsive acceleration phase reproduces the true curve in a 
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reasonable way. Higher time cadence data produce a better estimate of the acceleration 
duration. Somewhat surprisingly, the horizontal uncertainties on the acceleration peak are 
larger when a higher sampling rate is used (<2 min). This is due to the fact that small time 
steps between the HT data points with errors cause larger uncertainties in the subsequent 
derivative. Thus, the enhancement of the time cadence can only improve the acceleration 
profiles when also the errors in the CME HT measurements are reduced. 



4. CME acceleration compared to flare energy release 



The flare/CME events under study are: 2007 June 3 (C5.3), 2007 December 31 (C8.3), 
2008 March 25 (Ml. 7). Figures HI El and [10] show for the three events the evolution of the 
erupting CME in the low corona in EUVI 17lA and white light images from COR1. In both 
instruments we follow the leading edge of the CME as indicated by crosses in the figures. 
In EUVI 171 A (T ~ 1 x 10 6 K) typical CME features as the frontal rim and the cavity 



are observed (see also 



Aschwanden et al. 



20091 ). However, the embedded cooler prominence 



material (T ~ 1 x 10 4 K), visible in the subsequent coronagraph images, is missing. This 
is reasonable, since the contrast of promine nces against the c oronal background drops 



sharply for lines formed at T > 3 x 10 5 K (INoyes et al.l 



19721 ). The obvious similarity 



of the morphology and the temporal evolution of the developing/propagating CME 
str ucture justifies the 



sec 



Maricic et al. 



ink between EUV and coronagraph observations (for a discussion 
20041 ). By combining measurements of the leading edge of the erupting 
structure from EUVI and COR1 images we derive the kinematics of the CME with typical 
uncertainties of ±0.02-0.05 R for EUVI and ±0.12 R for COR1. We stress that for each 
event under study the location of the CME source (estimated by the flare position) is close 
to the limb (±10-25°). Thus, we expect that the influence of projection effects on the 



derived kinematics is small. We also note that for all events large-scale coronal waves were 
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detected in the EUVI 195A image sequences. 



4.1. 2007 June 3: C5.3 flare/CME event 

The 2007 June 3 CME event is ass ociated with a C5.3 GOE S class flare at heliographic 



position S08E67 (EUVI event catalog; 



Aschwanden et al 



2009). The position angle of 



STEREO-A with respect to Earth is 7°, hence for STEREO-A the CME was observed ~16° 
off the spacecraft plane-of-sky. The event was observed by EUVI with a high cadence of 
75 sec. 

Figure [5] shows the RHESSI HXR image (top) reconstr ucted in the energy band 



30-50 keV with the CLEAN algorithm f lHurford et al 



20021 ) together with its spectrum 



(bottom) integrated over 20 sec around the peak of the HXR emission. The spectrum was 
fitted with a thermal plus non-thermal model in the energy range 6-200 keV. The HXR 
source, located close to the eastern limb, is compact with some indication of two footpoints 
(but not clearly resolved). From the RHESSI spectral fit we derive a photon spectral index 
7 = 2.4 which indicates a very flat (i.e. hard) spectrum and a photon flux density at 50 keV 
during the event peak, F 50 = 0.93 photons s" 1 cm" 2 keV -1 . 

Figure |6] shows the CME distance, velocity, and acceleration time-profiles for the 
distance range 1.0-3.3 R Q (left panels from top to bottom). To enlarge the details on the 
impulsive acceleration phase, the right panels of Fig. O show the same curves but for the 
distance range up to 1.7 R . For comparison, the RHESSI HXR flux of the associated flare 
in the energy range 30-50 keV is overplotted in the bottom panels of Fig. |6j RHESSI reveals 
impulsive and powerful HXR emission with a distinct burst of short duration (less than 2 
min). The evolution of the acceleration profile of the CME and the evolution of the flare 
HXR flux are highly synchronized peaking at 09:27 UT with an uncertainty of ±1.25 min 
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(CME) and 09:27:14 UT (RHESSI 30-50 keV HXR flux). The CME reveals a large peak 
velocity (~1170 km s _1 ) and peak acceleration (~5.1 km s~ 2 ), which occur within the EUVI 
FoV, i.e. below 1.7 R . The CME accelerates within 8 min up to ~1170 km s _1 , then drops 
within the next 2 min down to a velocity of ~700 km s" 1 and further decelerates over the 
following 10 min until reaching a constant velocity of ~500 km s . If one would determine 
the CME speed solely from coronagraph observations, the velocity peak would be missed. 



4.2. 2007 December 31: C8.3 flare/CME event 



The 2007 December 31 CME event is associated with a C8.3 GO 



sourc e region lies behind the limb at roughly S09E102 ( IKrucker et al 



IS class flare. T 



lie 



2010 



Dai et al. 



20101 ) which makes this event special with respect to observations from Earth-view (GOES, 
RHESSI). From the RHESSI perspective the footpoints of the flare are occulted and we 
actually observe a coronal HXR source at non-thermal energies, with a good count statistic s 



in th e range 20-50 keV. This might be interpreted as Ma suda type source (IMasuda et al. 



19941 ) . For that day, MESSENGER flSolomon et al.ll200l[ ) is at a heliographic longitude 



of about E165 and observed the non-occulted flare in soft-X rays (similar to GOES). The 
flux measured by MESSENGER is about 2. 8 times higher than the flux observed by GOES 



which makes the flare actually of M2 class (IKrucker et al 



2010) 



Figure [H] shows the RHESSI HXR image reconstructed in the energy band 20-50 keV 
(top) together with its spectrum (bottom) integrated over 20 sec around the peak of 
the HXR emission and fitted in the energy range 6-90 keV. The image reveals that the 
strongest HXR emission comes from above the eastern limb, hence the footpoints of the 

?his m eans that the HXR flux entirely originates from the coronal source 



[e.g. 



are occulted. 


'his n 


Krucker et al. 


2008) 



20081 ). For this event we derive the RHESSI spectral fit parameters 



as 7 = 4.5, which means a softer HXR spectrum compared to the previous event, and 



13 



F^o = 0.17 photons s 1 cm 2 keV l . The considerably lower photon flux density and the 
softer spectrum as compared to the 2007 June 3 eve nt are not surprising since loop-top 



sources are generally less bright than the footpoints (IKrucker Sz Lin 



20081 ) 



At that day, STEREO-B had a position angle of 23° with respect to Earth, hence 
the CME was observed 11° off the STEREO-B plane-of-sky. Figure M shows the CME 
distance-time, velocity, and acceleration profile for the measured distance range 1-5 R (left 
panels). The closeup view (right panels) shows the distance range up to 2.5 R . We stress 
that for this event one data point (COR1) in the overlapping FoVs of EUVI and COR1 
(1.4-1.7 R Q ) reveals a shift between the both instruments, but which lies within the applied 
error bars. The CME reaches a maximum velocity of ~790 km s" 1 within <10 min, with 
a maximum acceleration of ~1.3 km s~ 2 . The impulsive acceleration phase of the CME is 
finished within the EUVI FoV and shows a peak at 00:50 UT (±2.5 min) which (within 
the uncertainties) is in accordance with the HXR evolution and peak time at 00:47:48 UT 
(Fig. [9j bottom panels). This is in particular interesting, as in this event we do not observe 
the HXR emission from the footpoints (which are occulted) but from a coronal source. 



4.3. 2008 March 25: Ml. 7 flare/CME event 



IS class flare. The 



20091 ) and the position 



The 2008 March 25 CME event is associate d with a M1.7 GO 
source region lies at S10E87 (EUVI event catalog; lAschwanden et al. 
angle of STEREO-B with respect to Earth is 24°, i.e. the CME propagates about 27° 
off the STEREO-B plane-of-sky. The event is observed with high cadence by EUVI and 
COR1 (75 sec and 5 min, respectively). RHESSI observed the peak phase of the flare HXR 
emission, but missed part of the rising phase due to spacecraft night. 



Figure [TT] shows the RHESSI HXR image in the energy band 18-30 keV (top) together 
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with its spectrum (bottom) integrated over 20 sec around the peak of the HXR emission. 
The spectrum was fitted in the energy range 6-100 keV, giving 7 = 3.5 and F 50 = 0.23 
photons s" 1 cm -2 keV -1 . The image shows a rather compact source of HXR emission 
revealing a weaker component in the corona in addition to the on-disk footpoint emission, 
suggesting that the second footpoint lies occulted behind the eastern limb (the energy 
range 18-30 keV is dominated by non-thermal emission, see spectrum in Fig. [TTj) . In 
EUVI 304A images a bright ejection is observed which may coincide with the coronal 
RHESSI source. 

Figure [TU] shows the sequence of STEREO-B EUVI running difference and COR1 
images. The images reveal a very thin and distinct spherically shaped front of the CME, 
reminding of the cross-secti on of a flux- rope torus. This shape is maintained in COR1 



white-light images (see also 



Aschwanden et al. 



20091 ) . The distance-time, velocity, and 



acceleration profiles are shown in Fig. [T2J overplotted with the RHESSI HXR flux in the 
energy range 18-30 keV. The COR1 data point lying in the overlapping FoVs of EUVI and 
COR1 fits well to the CME HT data measured by EUVI. Due to RHESSI night, the initial 
flare phase is missed, but the peak is clearly observed. RHESSI comes out of spacecraft 
night at ~18:44 UT revealing enhanced HXR emission with a maximum at 18:51:34 UT. 
The CME velocity starts to increase at 18:35 UT, the acceleration reaching a peak at 
18:50 UT (±3.5 min). We note that the acceleration time-profile shows two acceleration 
steps within the EUVI FoV. The CME reaches a peak velocity of 970 km s _1 within 
~20 min with an acceleration maximum of ~1 km s -2 . 



5. Summary & Discussion 

Both the spline fit and the regularization tool determine the acceleration phase of a 
CME (peak and duration) reasonably well. The main plus of the regularization tool over 
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the spline fit method is the provision of errors in time and amplitude. In all test cases, 
the true solution was included within this error estimation. We come to the somewhat 
surprising result that an increased image cadence (<2 min) does not reduce the uncertainty 
in determining the peak time of CME acceleration. This is due to the measurement errors 
of the CME leading edge. The inconsistencies of the CME leading edge as measured in 
different wavelength regimes (EUV and white light) as well as the change of intensity over 
the FoV of a single instrument do not allow us to narrow the measurement errors. The same 
holds for the acceleration amplitude for which the uncertainties lie in the range between 
10% and 50%. As a check for the goodness of the obtained velocity and acceleration profiles 
from the regularization tool, the normalized residuals between the regularized solutions and 
the distance-time measurements are investigated and found to be small (<1.5). 

Table [T] summarizes the relevant parameters derived from the CME and flare. From 
the three events under study we obtain that the peak acceleration of the CME is reached 
within a few minutes after its launch and at low distances of <0.4 R Q from the flare site, 
which we use as an estimate of the CME source region location. We would like to note that 
all events occurred close to the limb. Thus, projection effects are small, and the derived 
values should be close to the "true" ones. In these events, the impulsive acceleration phase 
is already finished before the CME is observed in the coronagraph COR1 FoV. The peak 
velocity of the CME is reached within heights of <2.2 R Q above the source region. The 
CME acceleration profile and the flare HXR emission are highly synchronized and peak 
almost simultaneously. The differences between the flare HXR peak and CME acceleration 
peak are in the range At < 2 min. Such differences lie within the limitations for deriving 
the CME acceleration peak time, which are caused by the finite cadence and measurement 
errors in the HT measurements. The uncertainties in the peak time of the CME acceleration 
estimated from the inversion method lie in the range 2-3.5 min (see Section [3]). 
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We note that there is no clear relation between CME velocity and GOES class of the 
associated flare. For the events under study these parameters vary only by a factor of 
~2. Statistical studies comparing the flare SXR peak flux and CME p eak velocity find a 



20021 ). 



linear correlation coefficient of about r ~0.47 which is not very strong ( iMoon et al. 
However, in our three events there seems to be a correlation indicated between the CME 
acceleration peak and the flare HXR peak flux, which both vary over a factor of 5 in the 
three events under study. Also some relation between the spectral slope of the HXR spectra 
and the CME acceleration is suggested (harder flare spectrum seems to be related to larger 
CME acceleration). Of course, three events are inconclusive for such a relationship but our 
results suggest that this should be tested on a larger event sample. 

Especially for the 2007 June 3 event intriguing results are derived. The single HXR 
burst of only ~2 min duration goes along with a rather high HXR flux density, CME peak 
acceleration (acceleration phase of ~8 min), and a very hard HXR spectrum. In contrast to 
that, a relatively low total energy in flare-accelerated electrons and low GOES classification 
is revealed. Further, the 2007 June 3 CME shows quite unusual behavior with respect to its 
very strong deceleration profile with values of —3 km s~ 2 and the velocity decreases from 
1100 km s" 1 to 500 km s" 1 very low in the corona ( <0.5 R^). Typical C ME deceleration 



Vrsnak et al 



2004aJ) and are due to 



values are in the order of —0.01 to —0.1 km s (e.g. 
the interaction with the solar wind flow. After the magnetic driving force of the eruption 
ceases, the CME is slo wed down by the drag force until its velocity adjus ts to the speed 



of the sol ar wind (e.g. 



therein) . 



Reeves 



Car gill 



2004 



Manoharan 



2006 



Vrsnak et al. 



2008 



and references 



( 120061 ) studied the relationship between the CME acc eleration and 



the 



thermal e nergy release rate in a loss-of-equilibrium flux rope model (cf. 



Lin &: Forbes 



2000 



Lin 



20021 ). They found that for high background magnetic fields and fast reconnection rates 



the evo l ution of the CME acceleration and flare energy release rate are well correlated. 



Reeves! ( 120061 ) also points out that due to rapid reconnection the formation of a long current 
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sheet in the wake of a CME is prevented and consequently the energy release is inhibited. 
For the 2007 June 3 event, we therefore assume a localized strong magnetic field from 
which the CME erupts as well as strong overlying fields. On the one side this yields high 
reconnection rates over a short time range, on the other side a strong overlying magnetic 
field would drastically decelerate the CME at low coronal heights. 

Well observed flare/CME events close to the limb, such as presented in this study, are 
less affected by projection effects when analyzing the CME parameters and enable us to 
reliably determine its kinematics and acceleration values. However, close to the limb the 
associated flare is observed either as "classical" on-disk HXR emitting source (footpoints 
of flare loops at chromospheric heights) or as coronal loop-top emission only, when the 
bright footpoint emission is occulted by the solar limb. In partially disk-occulte d events, 



loop-top sources could be studied on a statistical basis (cf. 



Krucker fc Lin 



20081 ) and it 



is most proba 



jle th at they are related to the initial location of particle acceleration (e.g. 



Krucker et al. 



20081 ) 



Battaglja fc Benzl fl2006f ) studied five RHESSI flares where HXR 



emission from both footpoints and the loop-top was observed. They found that the coronal 
and footpoint sources are well correlated with respect to their temporal and spectral 
evolution. Our results show in fact that the CME acceleration peak phase is well correlated 
with both types of HXR sources, i.e. coronal and footpoint HXR emission of the associated 
flare. 

In the present study we found for all three events under study, that the impulsive 
acceleration of the CME is finished at distances of <0.4 R Q above the solar surface, i.e. 
within the EUVI FoV and before reaching the COR1 FoV. Hence, from coronag raphic 



observ ation s only, the ear l y CM E evolution may be considerably misinterpreted. 



Vrsnak 



f l!990h and 



Chen fc Kralll ( 120031 ) obtained from analytical models, in which CMEs are 



treated as toroidal field structures, that the peak of the acceleration should be reached 
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when the half-separation of the CME footpoints is comparable to the CME height above 
the flare site. For the 2008 March 25 event the CME footpoint separation can be derived 
with adequate accuracy (see Fig. [7]). We obtained a CME half-separation of ~210±30 Mm, 
which is comparable to its peak height of ~280 Mm. 

Our results point to a possible relation between the flare HXR spectra and the CME 
acceleration (the event with the highest CME acceleration has the highest HXR flux density 
and flattest spectrum) whereas th e relation to the pea k power in electrons and the GOES 



soft X-ray peak flux seems weak. 



Maricic et al.l (120071 ) pointed out that the reconnection 



rate is more relevant for the CME acceleration than a strong heating and non-thermal 
particle acceleration. From our results it seems that the efficiency of accelerating particles 
to high energies is better correlated to the CME acceleration than the total number and 
energy in electrons (though we note that more statistics on this is needed). Such results 
may provide new constraints on electron acceleration models and its magnetic geometry as 
the particle acceleration mechanism has to go along with a rapid closing of magnetic field 
lines which is needed to drive the CME. 



6. Conclusions 



The present study based on high cadence data from STEREO-EUVI, STEREO-COR1, 
and RHESSI evidences a very close relation between flares and CMEs. From three well 
observed impulsive events that occurred close to the limb, we derive that the CME 
acceleration profile and energy release of the associated flare evolve in a synchronized 
manner. This supports the "standard" flare/CME model which predicts a feed-back 



relationship between the 
assoc iated flare (e.g 



arge-scale CME accelerat i on and the energy re 



the 


arge- 


Lin 


2004; 



Maricic et al. 



2007 



Temmer et al 



2008 



ease process in t he 



Vrsnak fc Oliver 



20081 ): After the magnetic structure looses equilibrium and starts rising, a current sheet is 



-In- 



formed below the rising structure (presumably a flux- rope), becoming a site of magnetic field 
reconnection. The reconnection has two important consequenc es for the CME acceleration: 



it reduces the downward-acting tension of the overlyin g field (ILin 



additional poloidal flux to t 
driver of the eruption (e.g . 



le flux rope 



Subramanian fc Vourlidas 



Chen 



1989 



( Vrsnak 


2008 


, which is 


Vrsnak 


1990 




Chen 


1996 



20041 ). and it supplies 



Kliem fc Torok 



2006 



20071 ). On the other hand, the upward moving CME drives mass 
inflow into the current sheet, driving further magnetic reconnection and particle acceleration 
(see the cartoon in Fig. [T3|) . Thus, reconnection directly relates the CME acceleration and 
the energy release in the associated flare, leading to the close synchronization of these two 
phenomena. 

The results from the presented study provide strong evidence for the feed-back 
mechanism between the flare energy release and the CME acceleration. In our study, all 
events were fast eruptions, hence, we can not draw conclusions on other types of eruptions 
(fa iled eru p tions or gradual CMEs without flares). According to model calculations 



by iReevesI (j2006l ) a good correlation (within 2 min) between the flux rope acceleration 
and thermal energy release rate is expected for fast reconnection events with high 
background magnetic fields. This should relate to impulsive and strong events, whereas 
no synchronization is expected for weak and gradual events. Finally, we note that high 
cadence observations of CMEs in non-coronagraphic images enable us to study in detail the 
CME impulsive acceleration phase revealing peak acceleration as high as ~5 km s -2 even 
for C-class flare events. We also note that the height of the CME peak acceleration in such 
impulsive events is much lower (<0.4 R Q ) than previously assumed. 
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Table 1: Summary of the CME and flare characteristics for each event. We give the date, the 
GOES flare class, the difference At [min] between CME peak acceleration and flare HXR 
peak, the derived CME peak velocity v max [km s" 1 ] and peak acceleration a max [km s~ 2 ] with 
errors (95% confidence level), the CME height h [R ] from the source region, the electron 
flux density at 50 keV, F 50 [photons s -1 cm" 2 keV" 1 ], the flare HXR spectral slope 7, and 
the total power in electrons P25 [10 27 erg s -1 ]. 

Date Class At v max (h) a max (h) F 50 7 P 25 

03-Jun-07 C5.3 0.1 1171±359 (0.48) 5.1±2.4 (0.26) 0.93 2.4 1.75 

31-Dec-07 C8/M2 2.0 785±125 (2.10) 1.3±0.4 (0.25) 0.17 4.5 4.09 

25-Mar-08 Ml. 7 1.5 967±173 (1.97) 1.0±0.2 (0.40) 0.23 3.5 1.88 
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Fig. 1. — Top: Synthetic curve of CME kinematics (black) together with the data (asterisks), 
results derived from the regularization tool (dashed blue) and the spline fit (pink). Middle: 
derived acceleration profile (black: true solution, blue: regularization tool result, pink: spline 
fit result). Bottom: normalized residuals (we note that in the left panel case the residuals 
make no sense due to the division by errors of the order of zero). Left: sampling rate is 
smaller than for the "true" curve, but the data are without noise. Right: same sampling 
rate as left but with noise added to the data. 
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Fig. 2. — Top to bottom: same as in Fig. [TJ Left: Impulsively accelerated CME evolution 
with rather high noise and errors. Right: Close up view for the time range during the peak 
acceleration phase. 
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Fig. 3. — Top to bottom: same as in Fig. [TJ Realistic scenario of CME evolution applying low 
noise to the data and error bars comparable to those from the observations. Left: Simulated 
is a time cadence of data points comparable to the actual observations. Right: same curve 
but with higher time cadence. 
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Fig. 4.— 2007 June 3, C5.3 flare/CME event observed with STEREO-A. Sequence of 
EUVI 171 A running difference images and COR1 total brightness images. Crosses indicate 
the measured leading edge from which the CME kinematics (shown in Fig. [6]) is derived. 
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Fig. 5.— 2007 June 3 C5.3 flare. Top: RHESSI 30-50 keV HXR image integrated over 20 sec 
around the flare peak using front detectors 2 to 8 (except 5 and 7) reconstructed with the 
CLEAN algorithm. Bottom: RHESSI spectrum and fit components for the energy range 6- 
200 keV. The isothermal fit is indicated as dashed-dotted line, the non-thermal power-law fit 
as dashed line, and the sum of both components as thick gray line. The derived parameters, 
electron flux density at 50 keV F^q [photons s _1 cm -2 keV -1 ], photon spectral index 7, and 
power P 2 5 [erg s _1 ] in electrons above a cutoff energy of 25 keV are given in the legend. 



-31 - 




Fig. 6.— CME kinematics for the 2007 June 3 C5.3 fiare/CME event. Top to bottom: CME 
distance, velocity, and acceleration against time together with the background subtracted 
flare HXR flux (red). Left: full height range, right: closeup view onto the early evolution 
phase. The pink dashed line in the right top panel (CME distance-time measurements) 
shows the regularized solution returned from the inversion technique. The grey shaded area 
in the CME velocity and acceleration curves indicates the 95% confidence level. 
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Fig. 7.— 2007 December 31, C8.3 (M2) flare/CME event observed with STEREO-B. Se- 
quence of EUVI 17lA running difference images and COR1 images. Crosses indicate the 
measured leading edge from which the CME kinematics is derived (Fig. [9]). 
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Fig. 8.— Same as in Fig. EI but for the 2007 December 31 C8.3 (M2) flare. The HXR image 
is integrated over the energy band 20-50 keV. 
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Fig. 9.— Same as Fig. M but for the 2007 December 31 C8.3 (M2) flare/CME event. The 
gray bar indicates the time range of RHESSI night (N). 
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Fig. 10.— 2008 March 25, Ml. 7 flare/CME event observed with STEREO-B. Sequence of 
EUVI 171 A running difference images and COR1 images. Crosses indicate the measured 
leading edge from which the CME kinematics is derived (Fig. fT2l) . 
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Fig. 11. — Same as in Fig. |5] but for the 2008 March 25 Ml. 7 flare. The HXR image is 
integrated over the energy band 18-30 keV. Contours show levels of 50%, 70%, and 90% of 
maximum emission for the energy band 30-50 keV. 
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Fig. 12.— Same as Fig. El but for the 2008 March 25 Ml. 7 flare/CME event. The gray bar 
indicates the time range of RHESSI night (N). 
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Fig. 13. — Cartoon illustrating the flare-CME feedback mechanism between the large-scale 
CME dynamics and small-scale flare processes. The upward moving CME evacuates the 
area in its wake, boosting mass inflow from aside into the reconnection region. The more 
mass and frozen-in magnetic field is transported into the region, the higher is the magnetic 
reconnection rate leading to larger flare energy release and to more efficient acceleration of 
particles. The successive closing of magnetic field lines due to reconnection increases the 
poloidal flux B^, in the eruption, which leads to a stronger upward oriented magnetic driving 
force (Lorentz force). 



